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Abstract. We study the existence and shape preserving properties of a generahzed 
Bernstein operator _B„ fixing a strictly positive function /g, and a second function fi such 
that /1//0 is strictly increasing, within the framework of extended Chebyshev spaces t/„. 
The first main result gives an inductive criterion for existence: suppose there exists a 
Bernstein operator Bn : C[a, 6] t/„ with strictly increasing nodes, fixing /o,/i £ Un- 
If Un C Un+i and Un+i has a non-negative Bernstein basis, then there exists a Bernstein 
operator Bn+i : C[a,6] Un+i with strictly increasing nodes, fixing /o and /i. In 
particular, if /o, /i, /n is a basis of Un such that the linear span of fo,..,fk is an 
extended Chebyshev space over [a,b] for each k = 0, ...,71, then there exists a Bernstein 
operator i?„ with increasing nodes fixing fo and /i. The second main result says that 
under the above assumptions the following inequalities hold 

Bnf > Bn+lf > f 

for all (/o, /i)-convex functions f G C [a,b] . Furthermore, Bnf is {fo, /i)-convex for all 
(/o, /i)-convex functions / G C [a, 6] . 



1. Introduction 

Given n E N, the space of polynomials generated by {1, x, . . . , x"} on [a,b] is basic 
in approximation theory and numerical analysis, so generalizations and modifications 
abound. However, from a numerical point of view it is a well known fact that the Bernstein 
bases functions Pn^k = {1 — x)^~^ behave much better and, in the sense of [2], provide 
optimal stability. The associated Bernstein operator _B„ : C [0, 1] — > ?7„, defined by 

k=o V / \ / 

has been the object of intensive research. As is well known, the polynomials Bnf converge 
to / uniformly although the convergence might be very slow. More important is the fact 
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that the Bernstein operator Bn reduces the variation and preserves the shape of /. In 
particular, if / is increasing then Bnf is increasing, while if / is convex then Bnf is convex, 
see e.g. [13J. And the derivative of Bnf of a function of class converges uniformly 
to /', cf. [22j, pg. 25. For this reason Bernstein bases and operators are fundamental 
notions. 

In Computer Aided Geometric Design (CADG) one is often interested, for instance, in 
rendering circumferences and other shapes not given by polynomial functions. It is thus 
natural to try to extend the preceding theory to more general spaces, containing not only 
l,x, . . . ,x"', but also, say, sine and cosine functions, while keeping as many of the good 
properties of Bernstein bases and operators as possible. If one generalizes the space of 
polynomials of degree at most n by retaining the bound on the number of zeros, one is led 
to the notion of an extended Chebyshev space (or system) f/„ of dimension n + 1 over the 
interval [a, b]: f/„ is an n + 1 dimensional subspace of {[a, b]) such that each f & Un has 
at most n zeros in [a,b], counting multiplicities, unless / vanishes identically. Recently, 
a rich mathematical literature has emerged concerning generalized Bernstein bases in the 
framework of extended Chebyshev spaces, see [9], [lO], [11], [12], [23], [21], [25], [26], [27], 

m, M, m- 

It is well-known that extended Chebyshev spaces possess non-negative Bernstein bases, 
i.e. collections of non-negative functions Pn,k, k = 0, n, in f/„, such that each p„ ^ has 
a zero of order A; at a and a zero of order n — k ai b, for k = 0, n. Assuming that Un 
has a non-negative Bernstein basis Pn,k, k = 0, n over the interval [a, b], it is natural to 
ask whether one may associate a Bernstein operator Bn '■ C [a, b] — > Un with properties 
analogous to the classical operator defined in ([1]). We consider operators Bn of the form 

n 

(2) Bn {f) = Y.f {t 

n.k ) Oin,kPn,k 

k=0 

where the nodes tnfi, ■■■,tn,n belong to the interval [a, b], and the weights anfi, are 
positive. But it is not obvious how the nodes and weights should be defined. Recall 
that the classical Bernstein operator reproduces the constant function 1 and the identity 
function x. We mimic this feature by requiring that Bn fix two functions /o, /i G Un, i.e. 
that 

(3) S„(/o) = /oand5„(/i) = /i, 

where throughout the paper it is assumed that /o > and that /i//o is strictly increas- 
ing, unless we explicitly state otherwise. Functions /o and /i satisfying the preceding 
conditions form a Haar system, in the terminology of [T^ pg. 25] . Following the Editor's 
suggestion, we shall call (/o, /i) a Haar pair. We show in Section 2 that after choosing 
/o and /i in Un, the requirements Bn (/o) = /o and Bn (/i) = /i, if they can be satisfied, 
uniquely determine the location of the nodes and the values of the coefficients; in other 
words, there is at most one Bernstein operator Bn of the form ([2]) satisfying ([3]) (observe 
that the only restriction on the nodes is that they belong to [a, b]). 
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The question of existence of a Bernstein operator in the above sense is studied in [T] 
and [2]. Here we present a new, inductive criterion for the existence of making this 
paper for the most part self-contained. Let /q, /„ G [a, b] and assume that for each 
k = 0, ...,n, the hnear space Uk := {fo,...,fk), generated by /o,...,/^, is an extended 
Chebyshev space of dimension k + 1. Then, for every k = 1, n, there exists a Bernstein 
operator Bk : C [a, b] — > Uk fixing /o and /i, whose sequence of nodes is strictly increasing 
and interlaces with the nodes of Bk-i, cf. Corollary [71 

Sections 3 and 4 deal with the shape preserving properties of the generalized Bernstein 
operator Bn- We shall utilize a generalized notion of convexity, (/o, /i)-convexity, which, 
according to [TB], p. 376, is originally due to Hopf, in 1926, and was later extensively de- 
veloped by Popoviciu, especially in the context of Chebyshev spaces. Ordinary convexity 
corresponds to (1, a;)-convexity. 

Assume there exists a Bernstein operator i?„ : C[a, b] —>■ Un fixing /o and /i. We shall 
show that if / G C [a, b] is (/o, /i)-convex, then 

Bnf > f, 

thus generalizing the same inequality for the standard polynomial Bernstein operator 
acting on convex functions. Assume next that Bn has strictly increasing nodes, that 
Un C Un+i, and that the latter space has a non-negative Bernstein basis. From the 
results in Section 2 we know that there exists a Bernstein operator Bn+i : C[a, b] — + Un+i 
fixing /o and /i. In Section 3 we show that 

Bnf > Bn+lf > f 

for all (/o, /i)-convex functions / G C [a,b] , generalizing once more the corresponding 
result for the standard polynomial Bernstein operator. In Section 4 we prove that under 
the preceding hypotheses, Bn preserves (/o, /i)-convexity, i.e., Bnf is (/o, /i) -convex for 
all (/o, /i)-convex functions f E C [a,b]. A similar result is obtained for the so-called /o- 
monotone functions /. These last results follow from the general theory of totally positive 
bases and their shape preserving properties. 

To put in perspective the inductive existence criterion indicated above. Section 5 (spe- 
cially. Theorem [23]) clarifies issues regarding the existence of "good" Bernstein operators, 
defined using non-decreasing nodes tn,k < ^n,fc+i in a suitable interval. We focus on the 
linear space generated by the functions cosx, and sinx on [0,6]. It is well known 

that normalized, totally positive bases (such as < (^)x'^ (1 — > in the polyno- 

mial case) posses optimality properties from the viewpoint of geometric design and shape 
preservation (cf. [6], [7], [T7]). 

It might be though that having such good bases would be enough to define a Bernstein 
operator fixing 1 and x. But this is not the case. The space U3 has a normalized, totally 
positive basis for every b G (0, 27r) (cf. [9], [23]). However, the existence of a Bernstein 
operator fixing 1 and x imposes the stronger condition b G (0,po], where po ~ 4.4934 is 



4 



J. M. ALDAZ, O. KOUNCHEV AND H. RENDER 



the first positive zero of 6 i-^ sin 6 — 6 cos 6: Wlien po < ^ < 27r it is not possible to find 
nodes in [0, b] so that 1 and x are fixed by the operator. When b G {tt,po], the nodes do 
belong to [0,6], but they fail to be non-decreasing. Thus, a Bernstein operator can be 
defined, but it lacks desirable properties; in particular, it does not preserve convexity. To 
ensure the existence of a "good", convexity preserving Bernstein operator, the stronger 
condition 6 < vr must be imposed, and to have a strictly increasing sequence of nodes we 
need even more: b < n. 

This paper is essentially self-contained. For simplicity, we consider only real valued 
functions when dealing with existence questions. Regarding shape preserving properties 
it is of course natural to consider real-valued rather than complex-valued functions. 

We thank the referees and the Editor, both for their very thorough reading of this 
paper and for their many suggestions, which lead to a substantial rewriting of the present 
article. 

2. Bernstein operators for Extended Chebyshev Spaces. 

We now introduce the concept of a Bernstein basis and of a non-negative Bernstein basis 
for a linear subspace f/„ C C" [a, b] of dimension n + 1. In the literature, the expressions 
"Bernstein like basis" or "B-basis" are often used instead of "Bernstein basis" . 

Definition 1. Let Un C C" [a,b] be a linear subspace of dimension n + 1. A Bernstein 
basis (resp. non-negative Bernstein basis) for Un is a sequence of functions (resp. non- 
negative functions) Pn,k, k = 0, ...,n, in Un, such that each Pn,k has a zero of exact order 
k at a and a zero of exact order n — k at b, for k = 0, n. 

By construction, (pn\{(^) ) is a triangular matrix with nonzero diagonal entries. 

V ' / i,k=0,...,n 

Hence, a Bernstein basis is indeed a basis of the linear space f/„. Furthermore, the basis 
functions are unique up to a non-zero factor, see e.g. Lemma 19 and Proposition 20 in 

As we indicated in the introduction, extended Chebyshev spaces always have non- 
negative Bernstein bases. To make this paper as self-contained as possible, we briefly 
indicate the reason: Let {Hq, be a basis for Un- To obtain a nonzero function pn,k 

with (at least) k zeros at a and (at least) n — k zeros at b, write pn,k '■= oo^o + ••• + cinhn- 
We impose the condition of having k zeros at a (which leads to k equations) and n — k 
zeros at b (which gives n — k additional equations). Having n + 1 variables at our disposal, 
there is always a non-trivial solution. The assumption that Un is an extended Chebyshev 
space guarantees that p„ ^ has no more than n zeros, so it has exactly k zeros at a and 
n — k zeros at b. In particular, pn,k is either strictly positive or strictly negative on (a, b). 
Multiplying by —1 if needed, we obtain a non-negative pn,k- 

In Proposition 3.2 in [26] it is shown that a subspace f/„ C [a, b] possesses a Bernstein 
basis Pn,k-, = 0, n if and only if every non-zero f & Un vanishes at most n times on the 
set {a,b} (and not on the interval [a, 6]). We mention that the existence of a Bernstein 
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basis in a space Un C [a, b] is a rather weak property; e.g. it does not imply the 
non- negativity of the basis functions Pn,k, k = 0, ...,n, nor the existence of Bernstein bases 
on subintervals [a,P] of [a,b]. 

The next two results are essential tools and standard techniques in CAGD in the context 
of degree elevation. 

Proposition 2. Assume that the linear subspaces Un C f/n+i C C^^^ [a, b] possess Bern- 
stein bases pn,k, k = 0, ...,n, and pn+i,k, k = 0, ...,n + 1. Then 

Pn+l,k l"J Pn+l,k+l 

for each k = 0, n. 



/ 4\ i^n,k\"'J ^n,k V"/ 
14J Pn,fc — —Pn+l,k + ;77Pn+l,fc+l 



Proof. Since G f/„+i, the function pn± is a linear combination of the basis functions 
Pn+i,k, k = 0, n + 1. Using the fact that pn,k has exactly k zeros at a and n — k zeros at 6, 

we see that Pn,k = aPn+i,k + PPn+i,k+i for some a, /3 G M. Then pJJ^]. = apjf+i + 
and inserting x = a yields 



Similarly, p^"^ = ap^J^^^l + jSp^^^^l^-^ and inserting x = b implies that 



Pn+l^l+l (b) 

□ 

Lemma 3. Under the hypotheses of the preceding proposition, assume additionally that 
the functions in the Bernstein bases are non-negative. Then 

(5) -Tpj > U and in-k') 

Pn+l,fc («) Pn+l,fc+l (^) 

/or eac/i = 0, n. 

Proof. If k = or k = n the assertion is obvious. If 1 < A; < n, then the first inequality 
in ([5]) can be obtained from (jl]): Divide both sides by Pn+i,k{x), and then let x I a. The 
second inequality follows in an analogous way. Alternatively, ([5]) can be derived, without 
using (jlj), from the well known and elementary fact that if / G C^'^^I) has a zero of order 
k at c, then 

(6) k\ ■ hm = /e^) (c) . 

^^'^ {x - c) 

Of course, the same formula holds for one side limits. □ 
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Let Un C [a, b] be a linear subspace of dimension n + 1 possessing a non-negative 
Bernstein basis Pn,k, k = 0, ...,n. We now introduce the concept of a Bernstein operator 
fixing a pair of functions: 

Definition 4. We say that a Bernstein operator Bn '■ C [a, b] — Un fixing the functions h 
and g exists, if there are points tnfl, ■■■,tn,n ^ [O', b] and coefficients anfl, > such 

that the operator Bn : C [a, b] —>■ Un defined by 

n 

(7) Bnf = ^f {tn,k) 0:n,kPn,k 

k=0 

has the property that 

(8) Bnh = h and Bng = g. 

We say that the sequence of nodes tn^, ■■■,tn,n G [a, is strictly increasing if 

tn,0 ^ tn^l ^ ... ^ tn,n- 

While the strict positivity of the coefficients an,k for k = 0,...,n is included in our 
definition of Bernstein operator, no restrictions are imposed on the nodes, save that they 
belong to [a, b]. For a natural example of a Bernstein operator without strictly increasing 
nodes, see Proposition [11] below, or Theorem [251 Section 5, in the case b = tt. However, 
if the nodes fail to be non-decreasing, then the Bernstein operator may behave in a 
pathological way, lacking convexity preserving properties (cf. Theorem [25]) . 

Two natural questions arise: when is the existence of a Bernstein operator guaranteed? 
and, is the Bernstein operator unique? It turns out that existence depends on additional 
properties of the space while uniqueness is easy to establish. 

In the development below we always assume that the Bernstein operator fixes a Haar 
pair (/o, /i) (i.e., /o G Un is strictly positive on [a, b] and /i G t/„ is such that the function 
/i//o is strictly increasing on [a, b]). The terms increasing and decreasing are understood 
in the non-strict sense. For a constant c positive means c > 0, while for a function / 
it means / > 0. Of course, once a Bernstein operator fixes a pair (/o, /i), it fixes every 
function in its linear span (/o, /i). 

We will consistently use the following notation. Assume that pj^k, k = 0, is a 
Bernstein basis of the space Uj. Given /o, /i G Uj, there exist coefficients Pj^, and 
■jjfi, such that 

j j 

(9) /o (x) = ^ f3j,kPj,k (x) and /i (x) = ^ lj,kPj,k (x) . 

k=0 k=0 

The next lemma answers the question of uniqueness positively. 

Lemma 5. Suppose that the linear subspace Un C C" [a, b], where n > 1, possesses a non- 
negative Bernstein basis {Pn,k}k=o- If there exists a Bernstein operator Bn : C[a,b] Un 
fixing the functions /o, /i G Un, then Pn,k > for all k = 0, . . . ,n, where the coefficients 
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Pn,k are given by Moreover, the nodes of Bn are defined, for k = and k = n, by 
tno = a and tnn = b, and in general, for k = 0, . . . ,n, by 



(10) t 



n,k • 



Jl\ / ln,k 



foj VA 



-'n,k 



with the '~fn,k given by Furthermore, the coefficients of Bn are defined, for k = 0, . . . ,n, 
by 

(11) «n,fc := yr: — 

JO[l^n,k) 

In particular, 

1 1 

(12) an,0 = r\ ^"^^ ^n,n = 77V 

Pn,o{a) Pn,n[0) 

Proof. Since 5„(/o) = /o and 

n 

(13) /o (x) = f3n,kPn,k {x) , 

fc=0 

we have 

n n 

/o (tn,k) OlkPn,k = Pn,kPn,k- 
k=0 k=0 

This entails that /o (tn,k) Oin,k = Pn,k, since {pn,k}k=o is a basis, and now ( ITTi) follows. 
Similarly, from -B„/i = /i and 

n 

(14) /i (2;) = ^-^^"'^ 

fc=0 

we obtain /i (tri,fc) Q^fc = 1n,k- Using /o > and q;„_a; > we see that Pn,k > 0. Dividing by 
/o (^n,fc) an,k = Pn,k, we find that t„,fc satisfies 

^j^g-j /l if'n,k) 1n,k 

fo (tn,k) Pn,k 

and now, since /i//o is injective, its inverse exists and we get ( ITOi) . Next, inserting x = a 
in (fT3!) and in (fT^ we obtain /o (a) = PoPn,o (a) and /i (a) = 7oPn,o (o). Thus 

fi (a) ^ 7n^ 
fo (a) /?n,0 ' 

and it follows by injectivity that t„.o = cl- An entirely analogous argument shows 
that tn,n = b. Since /o(a) = Pn,oPn,o{a) = fo{a)an,oPn,o{a) and /o(fe) = (3n,nPn,n{b) = 

fQ{b)an,nPn,n{b), foUoWS. □ 
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Lemma [5] tells us that to obtain a Bernstein operator Bn fixing /o and /i, the nodes 
tn^k must be the ones given by equation flTU]) . and the coefficients an,k by ffTTl) . A simple 
algebraic manipulation then shows that i?„ does fix /q and /i. To construct Bn, the 
difficulty lies in showing that for k = 0, ...,n, the numbers 

7n,fc 
Pn,k 

belong to the image of [a,b] under /i//o, so the nodes tn,k can be defined. Even if this 
is the case, it does not follow in general that the nodes are increasing (cf. Theorem [25]) . 
It does not seem trivial to characterize the spaces f/„ C [a, b] for which there exist a 
Bernstein operator fixing a predetermined Haar pair /o, /i G f/„, cf. [^. 

Here we present a new, inductive criterion: Existence of B^ with strictly increasing 
nodes entails existence of -Bn+i with strictly increasing nodes. Furthermore, the nodes at 
level n + 1 interlace strictly with the nodes at level n. 

Theorem 6. Suppose that the linear subspaces Un C Un+i C C"^^[a,6], where n > 
1, possess non-negative Bernstein bases Pn,k,k = 0,...,n, and pn+i,k,k = 0,...,n + 1 
respectively. If there exists a Bernstein operator Bn '■ C [a, b] —>■ Un fixing the functions 
/o,/i e Un, with strictly increasing nodes a = tnfl < < ••• < ^n,n = b, then there 
exists a Bernstein operator Bn+i : C [a, 6] Un+i fixing /o,/i, with strictly increasing 
and strictly interlacing nodes t„+i,o, ^n+i,n+i; that is, 

(16) a = tn+1,0 = tnfl < ^n.+l,l < ^n,l < < ^n,2 < " " " < ^n+l.n < ^n,n = ^n+l,n.+l = b. 

Proof Let us write /o = Y.ktl Pn+i,kPn+i,k and fi = Y^ktl 7n+i,fcPn+i,fc- By the preceding 
lemma, if the Bernstein operator Bn+i : C [a, b] — > f/n+i for (/o, /i) exists, then it has the 
form 

n+1 

Bn+lf '■= / (Wl,fc) Oin+l,kPn+l,k, 

k=0 

where the positive coefficients an+i,k are given by (fTTj) (with n + 1 replacing n), and the 
increasing nodes tn+i,fc are given by tn+i,o = o, by tn+i,n+i = b, and in general, by (ITOl) 
when A; = 0, . . . , n + 1. Thus, we need to show, first, that Pn+ifl, Pn+i,n+i > 0, in order 
to get the positivity of the coefficients an+i,k, and second, that 

(17) < 1^ < 1^ for A; = 1, ...,n, 

Hn,k-1 Pn+l,k Hn,k 

to obtain the (strict) interlacing property of nodes; note that 'jnfl/Pnfl = 1n+ifl/ Pn+ifl, 
since both quantities equal /i(a)//o(a), and similarly 7„,„//5„,„ = 7„+i,„+i//3n+i,n+i (since 
both quantities equal /i(6)//o(6)). 

At level n, by assumption the Bernstein operator is defined via the coefficients an,k > 0. 
Now the argument runs as follows: From the numbers an,k we obtain the f3n,k, and from 
these the Pn+i,k, which in turn give us the an+i,k- 
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Since f3n,k = fo{tn,k)an,k, it follows that (3n,k > 0. From /o(a) = Pn+i,oPn+i,o{a) and 
fo{b) = (3n+i,n+iPn+i,n+i{b) we See that Pn+ifi > and Pn+i,n+i > 0. We show next that 



for A; = 1, . . . , n, 



ns^ p _ f3 Pnji"') , ^ pi,k-i ^ jb) 

(18) fJn+l,k - Pn,k + fJn,k-l , 

from which the positivity of Pn+i,i, ■ ■ ■ , Pn+i,n follows by Lemma O Applying the index 
raising formula given by Proposition [2] to /o = J2k=o f^ri,kPn,k, we see that 



r _ a Pn,k (a) , 
Jo — Pn,0 T^Pn+1,0 + / ^ 

Pn+i,k (a) ^ 



Pn+l,k 



Pn,n{h) 
Pn+l,n+l [b) 

and we obtain (fT8|) . 

Regarding the interlacing property of nodes, another application of the index raising 
formula from Proposition [21 this time to /i = Yl^=oln,kPn,ki yields 

Pni{0) PnM-l^\b) 

(19) 7n+l,A: - 7n,fc (,) + 7n,/c-l 

for k = l,...,n. To show that < t^+ml^ equivalently, that 'ynk-i/3n+ik < 

Pn,k—1 Pn + l,fc ' ' 

7ri+i kPnk-i, we use formulas f|T8l) and f|T9l) to rewrite the latter inequality as 

(20) ' 

'"''" r"'%SM («) "^"''"ei-^^ (^) ; r'%s,. (a) ""'^"piti"^^ (^) ; "^""^ 

pC') (a) 

Simplifying and using -j^ > (by Lemma [3]) , inequality fl20l) is easily seen to be 

equivalent to < which is true by (ITOl) together with the assumptions that 

Pn , fc — 1 '7n , k 

/i//o is increasing and that tn,k-i < in,k- 

Inequality 1"^^''° < is proven in the same way. □ 

For the next corollary we do not a priori assume that /o > and /i//o is strictly 
increasing, since multiplying by —1 if needed, these properties can be obtained from the 
other assumptions. 

Corollary 7. Let /o,...,/„ G C"'[a,b], and assume that the linear spaces Uk generated 
by fo, fk are extended Chebyshev spaces of dimension k + 1 for k = 0, ...,n. Then for 
every k = 1, ...,n, there exists a Bernstein operator Bk : C [a,b] Uk fixing /o and fi, 
with strictly increasing nodes and strictly interlacing with those of B^-i- 
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Proof. Since Uq is an extended Chebyshev space over [a,b], the function /o has no zeros. 
Multiplying by —1 if needed, we may assume that /o > 0. Since Ui = (/o, /i) is an 
extended Chebyshev space over [a,b] it is easy to see that that /i//o is either strictly 
increasing or strictly decreasing. By multiplying /i by — 1 if needed one may assume that 
/i//o is strictly increasing. Let {pi,o;Pi,i} be a non-negative Bernstein basis for Ui. We 
define 

(21) BJ := ai_o/ (a) Pi,o + (b) pi^i, 

where ai^o = ^/pi,o{ci) and ai^i = l/pi^i{b). Since both functions (-Bi/o — /o) € Ui and 
(-Bi/i — /i) G t/i have a zero at a and another zero at b, and Ui is an extended Chebyshev 
space we see that these functions are zero, so Bi fixes /o and /i. And now the result follows 
by inductively applying Theorem [6] to each Uk+i in the chain Ui C U2 C ... C f/„. □ 

It is well known that given an extended Chebyshev space Un, one can find functions 
fo,---,fn € C"'[a,b] such that the linear spaces Uk generated by fo,...,fk are extended 
Chebyshev spaces of dimension k + 1 for k = 0, ...,n, see e.g. Proposition 2.8 of [26] (cf. 
also Definition 2.4 in [26]). The functions /o,---,/n can be constructed in the following 
way: first one shows that there exists a strictly larger interval [a, (3] D [a, b] such that 
Un is an extended Chebyshev system over [a, f3] (cf. [26l p. 351]). Take now n different 
points ^1, in the open interval (6, /?). For each k = 0, ...,n define a non-zero function 
fk e Un vanishing on ^i,...,^n-k- Then the linear spaces Uk generated by fQ,...,fk are 
extended Chebyshev spaces over [a, b] . The disadvantage of this procedure is that the 
choice of the functions /o and /i cannot be specified in advance, but does depend on the 
space Un- 

Thus, Corollary [7] implies the next result: 

Corollary 8. Let n > 1 and let Un be an extended Chebyshev space over [a, b] . Then it is 
possible to find a Haar pair /o, /i G f/„ and a Bernstein operator Bn : C [a, b] — > Un with 
strictly increasing nodes, such that Bn fixes fo and fi . 

Remark 9. Observe that the hypothesis of Corollary [S] is weaker than that of Corollary 
[TJ and so is the conclusion, since /o and fi are chosen a posteriori, cf. also the discussion 
at the beginning of Section 5. 

Specializing the preceding results to the case of exponential polynomials, the conclusions 
we obtain in the real case are stronger than those from [1] (however, [1] deals with the more 
general complex case). The space -E'(Ao,...,a„) of exponential polynomials with exponents 
Ao, . . . , A„ G C is defined by 

Exponential polynomials provide natural generalizations of the classical, trigonometric, 
and hyperbolic polynomials (see [33]), and the D-polynomials considered in [15] and [29] . 
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They also furnish (under additional assumptions) examples of extended Chebyshev spaces 
recently studied in CAGD for the purpose of representing parametric curves, see (TU], [12], 
[23]. Note in particular that when = for < A; < ri, A„+i = i and A„_|_i = —i, we 
have -E'(Ao,...,A„+2) = (1; • • • , 2;", sinx, cosx), where {gi, . . . ,gk) denotes the vector space 
spanned hj gi,...,gk. 

It is well known that -E(Ao,...,a„) is an extended Chebyshev space over any compact 
interval [a, b], for any choice of Aq, A„ G M (this is not true if some of the exponents are 
complex, but here, we only consider real eigenvalues A^; the reader interested in complex 
exponents may want to consult [IJ, [2] and [21] )• 

The results of the present paper give a simple proof of the existence a Bernstein operator 
with strictly increasing and strictly interlacing nodes, fixing a suitable Haar pair. Let us 
emphasize that the interlacing property of nodes does not follow from [1]. 

Theorem 10. Let Aq, A„ be real numbers, let a < b, and for k = 1, . . . ,n, let Uk '■= 

-E'{Ao,...,Afc) over [a,b]. If Xq = Xi, we set fo{x) := e^°^ and /i(x) := xe'^°^ , while if Xq 7^ Ai, 
we set fo{x) := e^^'^^^^O'^i^ and fi{x) := e''^'"^^^^'^^^ Then for each k = 1, ...,n, there is 
a Bernstein operator : C [a, 6] — >• E(^Xo,...,Xk) /j^^^^fl' /o one? fi, such that its sequence of 
nodes is strictly increasing, and those at level k > 1 strictly interlace with the nodes of 
Bk-i- 

Proof. It immediately follows from (1221) that -E(Ao,...,a„) is invariant under permutations of 
the eigenvalues Aj, so we can assume, with no loss of generality, that Aq < Ai. Note that 
if Ao = Ai, then Ui = E^XoM) = (/o,/i)- Whether we have Aq = Ai or Aq < Ai, /o > 
and /1//0 is strictly increasing. Since the spaces Uk = -E(Ao,...,Afc) over [a,b] are extended 
Chebyshev spaces for each k = 0,...,n, the result immediately follows from Corollary 

El □ 

We finish this section with a proposition illustrating our methods in the classical poly- 
nomial case. The article [20] exhibits a sequence of positive linear operators converging 
to the identity on C[0, 1] and fixing 1 and x^. This sequence is obtained by replacing x in 
([1]) with a suitably chosen function r„(a;) such that lim„rn(a;) = x. It is also possible to 
fix 1 and by using the generalized Bernstein operators considered here. As a matter 
of fact, it is possible to fix fo{x) = 1 and fi{x) = x^ for any j > 1 we wish (of course, if 
j = 1 we have the standard case). From Lemma [5] we know how to determine the nodes 
and the coefficients, i.e., how i?„ must be constructed. 

On the other hand, we cannot use CorollaryElto conclude that such a Bernstein operator 
Bnflj exists (the subscripts and j refer to the exponents of the functions being fixed) 
since whenever j > 1, the space Ui = {l,x^) is not an extended Chebyshev space over 
the closed interval [0, 1]: x^ has a zero of order j. And unlike the situation considered in 
Theorem [6], the sequence of nodes we obtain is not strictly increasing: given 1 < J < n, 
it is easy to see that t„,o = ■ ■ ■ = tn,j-i = 0, simply by counting zeros at a = 0, or by the 
argument given below. 
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Proposition 11. Fix j > 1, and let Un he the space of polynomials over [0,1] of degree 
at most n. For every n > j , there exists a Bernstein operator Bn^j '■ C[0, 1] Un that 
fixes 1 and x\ and converges in the strong operator topology to the identity, as n ^ oo. 
The operator Bn^j is explicitly given by 

^^■».^/'->-|:/(( :!::l!:::|::^:i! )'l(:)-''^--)- 

Proof. For the purposes of this argument we set Pn,k{x) := (^)x^(l — (this differs 

from the notation used in the introduction for the classical Bernstein polynomials, but 
it is more convenient here). The condition 1 = Bnfljl{x) = Ylk=o ^ri,kPn,k entails that 
= 1 for all n, k. We use the equality x^ = Bn^jx^ to determine the nodes tn^k- 
Writing 

n 

2'"' ^ ^ '~^n,kPn,ki 
k=0 

by (ITUI) we have t„_fc = 7^^. The coefficients 7^^^ can be obtained in several ways. A 
rather direct one follows next: 



= X- 

k=0 



k 



V /^^ - A ^fcf 1 _ ^)n-k _ ^ k{k-l)---ik-j + l) f n\ , _ 
eU-j7 ^ ^ -£-n(n-l)---(n-j + l)W ^ ^ 



k=j ^ ^ k=j 

Thus, 'jn,k = if < /c < j and 

kik-l)---{k-J + l] 



n{n — 1) • • • (n — j + 1) 

when j < k < n. Originally, we found 7„,fc using Lemma [5l but an anonymous referee 
tells us that equation ( l23l) is well known and can be obtained by blossoming (cf. [32] 
for an introduction to blossoms). And the Editor adds that the coefficients can also be 
determined via the dual functionals for B-splines. Apparently, though, the coordinates 
'jn,k of with respect to the Bernstein basis had not previously been used as we do here, 
to define a Bernstein operator fixing 1 and xK 

Next we prove convergence. For / = 1, . . . , j — 1, the inequalities 

- ? + 1 k-l k 
< < - 

n n — I n 

can be checked by simplifying and inspection. It follows that ((fc — j + l)/^)"* < tnk ~ 
ln,k < {k/ny , or equivalently, that < k/n — tn,k < (j — Thus, Bn^jx"" converges 

uniformly to x"^ for m = 0, 1,2, and by Korovkin's Theorem, Bnfljf f uniformly for 
all/eC[0,l]. " □ 
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3. Generalized convexity 



Let Bn denote the classical Bernstein operator defined in ([T]). W.B. Temple showed in 
3^ that for a convex function / the following monotonicity property 



(24) 



XI 



holds for all x G [0, 1] . In ^4j O. Arama proved that 



BJ (x) - (x) 



x(l 



N n— 1 

-T 

n(n+l)^ 



k k+1 k+1 
n' n + 1 ' n 



n 



k 



x'^il 



\ n— 1— fc 



where ^^^j / is the divided difference of second order, thus providing a simple 

proof of Temple's result. A similar formula (see Theorem 7.5 in [19]) is due to Aver- 
bach. We obtain analogous results for the generalized Bernstein operators considered 
here. These generalized Bernstein operators -B„ fix /q and /i instead of 1 and x, so rather 
than (1, x) -convexity, which is equivalent to standard convexity, the adequate notion for 
our purposes is (/o, /i)-convexity, to be defined next. We shall see that for (/o, /i)-convex 
functions /, the following holds: Bnf > -Bn+i/ > /• 

Next we specialize the definition given in ^31 p. 280] for an arbitrary number n of 
functions to the case n = 2. Later we will specialize it to the case n = 1. 

Definition 12. Let E cM.. A function f : E —>■ is called convex on E with respect to 
a Haar pair (/o, /i) if for all xq, xi, X2 in E with xq < xi < X2, the determinant 



(25) 



Det 



X0,Xl,X2 



fo (Xq) /o (Xi) /o (X2) 

(/):=det| /i(xo) /i(xi) /i(x2) 
/ (a^o) / (xi) f (xa) 



is non-negative. We shall also use the shorter expression "{fo, fi)- convex" . Likewise, we 
say that f is {fo, fi) -concave if -f is {fo, fi) -convex, and {fo, fi)-affine if f e (/o,/i). 



Remark 13. Note that the condition Det 



X0,Xl,X2 



(/) > for all xo,xi,X2 in E with 
Xq < Xi < X2 is equivalent to the same requirement but with Xq < Xi < X2. Of course, in 
the degenerate case Xj = Xj+i the determinant is zero, so it makes no difference whether 
or not this possibility is included in the definition. In other words, only the ordering of 
the points Xo,Xi,X2 actually matters. 

One of the standard definitions of convexity stipulates that the graph of / must lie 
below the segment joining any two given points on the graph. It is well known that an 
analogous characterization holds for (/o, /i)-convex functions, but with affine functions 
being replaced by {fo, /i)-affine functions. More precisely. 

Proposition 14. Let {fo, fi) be a Haar pair. Denote by ipl^ the unique function in 
Ui := {fo, fi) that interpolates f at the points xq < X2, xo,X2 G [a,b], i.e., '?/'4,a;2 (^0) ~ 
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/ (xo) and '?/'4,a;2 (^2) = / {X2) ■ Then f is (/o, f\)-convex if and only if for all xo,x,X2 
such that a < Xq < X < X2 < b, 

(26) / (x) < ^i^^^ (x) , 
and in this case, for all y G [a, b] \ [xq, X2], 

(27) / (y) > V'i,.. (y) ■ 
Proof. Let Detxo,x,x2 (/) be as defined in (l25l) . Observe that 

^^^Xo,X,X2 if) = ^^^Xo,X,X2 (/ ~ ^^^0,2:2) 

= - (/ (x) - V^4,^2 (^)) (/i (^2) /o (a;o) - /o (3^2) /i (xo)) • 
Since /1//0 is strictly increasing, /i (^2) /o (a;o) - /o (2^2) /i (a^o) > 0, so Det^o,a;,x2 (/) > 
is equivalent to / (x) < ipl^ X2 (•^)- 

Next, assume that (HEh holds for all Xo,x,X2 such that a < xq < x < X2 < b. Suppose 
that for some u G [a, 6] \ [xo,X2] we have f{u) < ^"4x2 i'^)- Without loss of generality 
we may assume that X2 < u. We interpolate between xq and u to obtain a contradiction: 
'^Lx2 (^0) = iplo,u i^o), while V4,X2 (^) > ^^4,^ (^)- Since ipf^^^^ - V'4,« has exactly one 
zero, it follows that ?/'4^^^ > ?/'4^„ on (xq, b]. In particular, /(xa) = ^/'4,X2 (^^2) > V'4,u (^2), 
which is impossible by f l26l) applied to / and V'4,tf '-' 

Note that by inequality (!26l) . convexity is the same as (1, x)-convexity. The next result 
generalizes to (/o, /i)-convex functions, the familiar inequality i?„/ > / for the classical 
Bernstein operator acting on convex functions. Here it is not assumed that -B„ is defined 
via an increasing sequence of nodes; it is enough to know that t„_fc G [a, b]. 

Theorem 15. Assume that for some n > 1, there is a Bernstein operator Bn fixing /o 
and fi. Then for every (/o, fi) -convex function f G C[0, 1] we have Bnf > f ■ 

Proof. Suppose B^ exists for some n > 1, and let £ > 0. We show that for an arbitrary 
X G [a, 6], Bnf{x) > f{x) — e. Assume that x G (a, 6) (the cases x = a and x = b 
can be proven via obvious changes in the notation, or just by using continuity). First, 
select 6 > such that BnS < e. Next, by continuity of /, choose h > so small that 
[x - h,x + h] C [a, b] and ■i/'f,^ < f + 6 on [x - h, x + h]. Then V'f./i x+h < / + ^ on 
[a,6] by (127D, so 

5„/(x) > 5„ (iji_h^x+h - ^) (^) = Bniji_h^x+h{^)-BJ{x) > ^f_^^^+^(x) -£ > f{x)-e, 
where for the last inequahty we have used fl26|) . □ 

We shall use below the following characterization of (/o, /i)-convexity, due to M. Bessen- 
yei and Z. Pales (cf. Theorem 5, p. 388 of [5J). While the result is stated there for open 
intervals, it also holds for compact intervals. Note that what we call here a Haar pair 
(i.e., /o is strictly positive and /1//0 strictly increasing) is called in [5] a positive regular 
pair. 
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Theorem 16. Let (/o, /i) he a Haar pair and let I := (/i//o)([o., b])- Then f G C[a, b] is 
(/o, fi)-convex if and only if (///o) o (/1//0) ""^ ^ C (/) is convex in the standard sense. 

Example 17. Consider the Bernstein operator -Bn,o,i fro^i Proposition [H], defined on 
C[0, 1] and fixing 1 and . It is easy to see from Theorem [16] that for s G (0,j), the 
function x^ is (1, x-^) -concave, while if s G (j, 00), x^ is (1, X'')-convex. Therefore, by 
Theorem [T5| for all x G [0, 1] we have BnfijX'^ < x'^ if s G (0, j) and Bn^jX^ > x^ when 
s G (j, 00). 

Our next objective is to obtain an analog of Arama's result (presented at the beginning 
of this section) for generalized Bernstein operators Bn- Here the interlacing property of 
nodes is used in an essential way. 

Proposition 18. Let Sk, Sk+i, Sk+2 £ [a-, &] be such that Sk < Sk+i < Sk+2, o-nd assume 
that Gk : C [c, d] ^ W is a functional of the form 

Gk (/) = akf (sfc) + bkf (sfe+i) + Ckf (sfc+2) , 

satisfying Gkifo) = Gk{fi) = 0. Then bk > if an only if Gk{f) < for all (/o, 
fi) -convex functions f & G [tk,tk+2]- 

Proof. Let / be (/o, /i)-convex, and let V'4,Sfc+2 function in (/o, /i) that interpolates 

/ at the points Sk and Sk+2- By f l26l) . 

if) =Gk{f- ^PU^,) =bk{f- (^^+1) < 

if and only if bk > 0. □ 

Theorem 19. Under the same hypotheses and with the same notation as in Theorem\^ 
let the linear functionals Gk., k = 1, ...,n — 1, be defined by 

rl-^'^ in) J.n+1-k) /,x 

^k [J) — J [l^n,k) 0:n,k—^i — — / [tn+l,k) Oln+l,k + J \tn,k~l) ^n.fc-l (n+l-k) 

Pn+l,k W Pn+l,k {bj 

Then 

n 

Bnf - Bn+lf = ^ Gk (/) ■ Pn+l,k- 

k=l 

In particular, if f is (/o, fi)-convex then B^f — -Bn+i/ > 0. 
Proof. Recall that 

n n+1 
n.k ) Oin,kPn,k 

and Bn+lf = ^ / (Wl,fc) <^n+l,kPn+l,k 
k=0 k=0 
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where tn,o — tn+i,o — a, tn,n — tn+i,n+i — b, and tn,fc-i < tn+i,k < tn,k for k — l,...,n. 
Using Proposition [2] we obtain 

p^''^ (a) 

n.k) (^n,k ^^-j Pn+l,k 
k=0 Pn+l,k W 

Oin,k (J'l^) — Pn+l,fc+l — E / ttri+l,fcPri+l,fc- 

fc=0 Pn+l,k+l (") fc=0 

It foUows from (fT2l) that the first summands (corresponding to = 0) of the first and the 
last sum are the same, so they cancel out. Likewise, the n-th summand of the second sum 
and the (n + l)-st summand of the last sum cancel out. Thus 

Bnf — Bn+if = 



fc=l 



r /, X Pn\ (^) r u N , f /J. \ Pn,k-l ^ (^) 

J V-n^k) Oln,k—^ — - / \tn+l,k) On+l.fc + / \tn,k-l) "n.fc-l (n+l-k) /,n 

^'n+l.fc (^) Pn+1,A: (^) 



Finally, let / be (/o, /i)-convex. Taking Sk = tn,k-i,Sk+i = tn+i,k, and Sk+2 = tn,k in 
Proposition [m we get Bnf — -Bn+i/ > 0. □ 

A very natural question, not touched upon here, is under which conditions a sequence 
of Bernstein operators for (/o, /i) converges to the identity. It follows from Theorems [15] 
and [I9] that if / is (/o, /i)-convex, then the sequence {Bnf}'^^i monotonically converges 
to some function g > f (assuming that a sequence of functions /o, /i, /2, ... are given such 
that (/o, /n) is an extended Chebyshev space of dimension n + 1 for each G N). But 
we have not determined which conditions will ensure that g = f ■ In this regard, we expect 
the strict interlacing property of nodes to be useful, since it entails, in a qualitative sense, 
that the sampling of functions is not "too biased". 

4. Total positivity and generalized convexity 

Let En : C [a,h] ^ f/„ be a Bernstein operator for the pair (/o, /i). In Section 3 we 
proved that Bnf > f for all (/o, /i)-convex functions / G C [a, 6] . This did not require an 
increasing sequence of nodes; it was enough to know that t„ G [a, h]. 

In this section we show that Bnf is (/o, /i)-convex for every (/o, /i)-convex function 
/ G C[a,6], provided that the nodes t„0)---)^n,n are increasing and Un is an extended 
Chebyshev space over [a, h] , and a similar result holds for the so-called (/-monotone func- 
tions. These statements follow directly from the more general results presented in [T9] 
concerning shape preserving properties of linear transformations with totally positive ker- 
nels. The connection between total positivity and shape preserving properties of bases is 
a classical subject that has been widely described, see e.g. [19] or the more recent survey 

n 
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The following definitions come from [19]. Let X and Y be subsets of M. A function 
K : X X Y M. is called sign- consistent of order m if there exists an 6^ G { — 1, 1} such 
that 



(28) Emdet 



( K{xi,yi) K{xi,y2) ... K{xi,ym) \ 
K{x2,yi) K{x2,y2) ■■■ K{x2,ym) 

\K{Xm,yi) ••• •••• K{Xm,ym) J 



> 



for all xi < X2 < ... < Xm in X and yi < y2 < ■■■ < ym in Y. If = 1 we shall call 
K positive of order m. A function K is totally positive if it is positive of all orders m 
with m G N, m > 1. Similarly, if one has strict positivity in fl28|) then K is called strictly 
sign- consistent of order m, and if in addition = 1, then K is strictly positive of order 
m. Strict total positivity means that K is strictly positive of all orders m G N, m > 1. 

The following result is well-known, see e.g. [26l p. 358], or the proof presented in [121 
pp. 342-344]: 

Theorem 20. Let Un C [a, b] be an extended Chebyshev space over [a, b] and let 
Pn,k, k = 0, ...,n, be a non-negative Bernstein basis for [a, b] . Then K : [a, b] x {0, n} 
M, defined by 

(29) K (x, k) := pn,k (x) , 

is totally positive, and K is strictly totally positive on {a,b) x {0, ....,n} . 

Following the notation of [9], [10], [23], we can deduce from the previous result that 
a non-negative Bernstein basis of an extended Chebyshev system over [a, b] is totally 
positive on [a,b], i.e., a B-basis. 

We cite from [T^l p. 284] the following result (specialized to the case of two functions 
Fq, Fi instead of a family Fi, F^). 

Theorem 21. Let X and Y be subsets o/ M, let Fq, Fi be functions on Y and let K : 

X X y — i> M 5e continuous, and positive of order 3. Let ^ be a non-negative sigma-finite 
measure and '■ C (Y) C (X) be defined by 



Bk{F) (x) := j^K{x,y)F{y)dii{y). 



If F is {Fq, Fi)- convex, then Bk{F) is {B^Fq, BKFi)-convex. 

From this we conclude: 

Theorem 22. Let Un be an extended Chebyshev space over [a,b] . Assume there exists a 
Bernstein operator Bn '■ C [a, b] — > Un fixing /o and fi, with increasing nodes tn,o < •••• < 
tn,n- If f ^ C[a,b] IS (fo, fi) -convex, then 5„ (/) is (fo, fi)- convex 
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Proof. Put X = [a, b] and Y := {0, n}. Define the function ip : Y ^ X hj (f (k) := tn,k, 
for k = 0,...,n. Observe tliat if is monotone increasing (thougfi perfiaps not strictly), 
so it is order preserving. Next, set /i := J22=o^n,kSk, where the an,k are the positive 
coefficients defining Bn and Sk is the Dirac measure at the point k G {0, ...,n}. With 
K{x, k) := Pn,k{x), we obtain, for every F & C (Y), 

» n 

Bk{F){x) := / K{x,y)F{y)dii{y) = y^F{k)an,kPn,k{x). 
■'y k=o 
Now let / G C (X) , and define F := f o ^ e C {Y) . Then 

n 

(30) Bk if 0(^){x) = ^f {tn,k) 0;n,kPn,k (x) = 5„ (/) (x) . 

A:=0 

If / G C (X) is (/o, /i)-convex, then F = f o ip is {f^ o ip, fi o ip)-convex, since ip preserves 
order (cf. Remark [T^ . Putting Fj = fj o (p for j = 0, 1, an application of Theorem fill 
shows that Bk (F) is [BkFq, BKFi)-coiivex. By formula fl30l) and the property that Bn 
fixes /o and /i one obtains 

BrFj = Bn ifj) = fj 

for J = 0, 1. Thus Bk (F) = Bn (/) is (/o, /i)-convex. □ 

In a similar way it is possible to obtain generalized monotonicity properties of the 
Bernstein operator. 

Definition 23. Let g > 0. We say that f is g-increasing on [a,b] if f / g is increasing on 
[a, 6] , I.e. if 



det 



g{xo) g{xi) 

for all xo < xi in [a, b] . 

The notions of (/-decreasing and (/-monotone are the obvious ones. Now / : [a, 6] ^ S 
is (/-increasing on [a, b] if and only if for all xq, Xi in [a, b] with xq < xi the determinant 

g{xo) g{xi) 
fixo) f{xi] 

is non- negative. But this condition is just the definition of convexity with respect to one 
function (the function g), obtained by specializing to = 1 the definition given in [T9l p. 
280]. Specializing also Theorem 3.3 in [19, p. 284] to one function (cf. Theorem [2T] above 
in the case of two functions) and arguing as in the proof of Theorem [22l we obtain: 

Theorem 24. Let Un be an extended Chebyshev space over [a, b] . Assume that there exists 
a Bernstein operator Bn '■ C [a, 6] Un with increasing nodes tnfl ^ •••• ^ ^n.n? fixing /q 
and fi. If g & (/o^/i) is strictly positive and f G C [a, 6] is g-monotone, then Bn (/) is 
g-monotone. 
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One of the referees points out that Theorems [24] and [22] can be proven without referring 
to Theorem 3.3 in p. 284], using instead elementary shape preserving properties of 
totally positive bases, as described in the surveys [8] and [17] . Moreover, Theorems [2l] 
and [22] in the special case /o = 1 and /i (x) = x follow immediately from [T71 Corollaries 
3.7 and 3.8, p. 162]. 

5. Normalized Bernstein bases and existence of Bernstein operators 

Let bn,k, = 0, n, be a basis of a given subspace Un d C [a, b] of dimension n + 1. The 
basis bn^k, k = 0, n is totally positive if the kernel K [x, k) := bn^k [x) is totally positive 
(in particular, the functions bn,k are non-negative). Suppose now that Un contains the 
constant function 1. Then the basis bn,k, k = 0, ...,n is called normalized if 



for all X e [a, b] . Normalized totally positive bases are important in geometric design due 
to their good shape preserving properties. J.-M. Carnicer and J.-M. Pena have shown 
that normalized, totally positive Bernstein bases are optimal from the shape preservation 
viewpoint, see [6], [7], [17]. Moreover it was proven in [9], and independently in [26], that 
a subspace Un of C" [a, b] of dimension n + 1 containing the constant functions, possesses 
a normalized, totally positive Bernstein basis, provided that both Un and the space of all 
derivatives Un '■= {/' : / € Un} are extended Chebyshev spaces over [a, b] . 

From this point of view it is natural to conjecture that to define a well-behaved Bern- 
stein operator (with increasing nodes) fixing /q and /i, it is enough to assume that Un 
possesses a normalized totally positive Bernstein basis and that (/o) , (/o, /i) are extended 
Chebyshev systems. However, we shall show by a counterexample that this is not true. 
We refer to [2] for a more detailed discussion under which conditions there might exist a 
Bernstein operator fixing a pair /o, /i G t/„. 

Consider the linear space U3 := (1, x, cosx, sin x) over the interval [0,6]. This space 
has been previously studied by several authors, see the references in [9] or [23]. It is 
well known that U3 and W3 (the space of all derivatives of functions in U3) are extended 
Chebyshev spaces over [0, b] for every b G (0, 2n) . Thus U3 possesses a normalized totally 
positive Bernstein basis for every b G (0, 27r) . By [H] this entails that the critical length 
of W3 for design purposes is 27r. However, we show in Theorem [2S] that for b sufficiently 
close to 271 (say, b > 4.5) there is no Bernstein operator from C[0, b] to W3 fixing 1 and x. 

The obstruction for employing Corollary [7] is due to the fact that neither (l,x, cosx) 
nor (l,x, sinx) are extended Chebyshev spaces over [0,6] for all 6 < 27r (for instance, 
sin X — X has a zero of order 3 at 0) so the chain of nested spaces cannot be continued 
beyond Ui = {l,x). By Corollary [8] it is nevertheless possible to construct a Bernstein 
operator fixing some Haar pair of functions go,gi G W3 = (1, cosx, sinx). Hence, by 



n 




k=0 
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Theorem [6] there is a corresponding Bernstein operator from C[0,6] to W3, fixing go and 
Qi, with strictly interlacing nodes. 

Theorem 25. Given b G (0,27r), let po be the first positive root of b ^ sinb — bcosb, 
(po ~ 4.4934/ Let W3 = (1, a;, cosx, sinx), /o = 1, and fi{x) = x. Then for every 
b G (0,po] there exists a Bernstein operator B3 : C[0,b] —>■ fixing 1 and x. The nodes 
to {b) , ti (6) , t2 {b), and t^ (6) satisfy the following inequalities: 



The operator B3 preserves convex functions whenever b G (0, tt] , but not for b G (vr, po] • 
Finally, if b & (po,27r), then there does not exist a Bernstein operator fixing 1 and x. 

To keep computations simple we shall present first two general propositions and a 
definition: We say that a subspace f/„ C C [a, b] is symmetric if / G implies that the 
function F defined hj F (x) := f {a + b — x) is in f/„. 

Assume that Un C C" [a, b] is a symmetric, extended Chebyshev space over [a, b] . Let 
Pn,k, k = 0, ....,n be a non-negative Bernstein basis of f/„, and let Pq, ...,Pn and 70, ...,7n 
be constants such that 1 = J2k=o PkPn,k (x) and x = Ylk=o 1kPn,k (x) for all x G [a, b] . 

Proposition 26. Suppose that Un C C" [a, b] is a symmetric, extended Chebyshev space 
over [a, b] containing the constant function /o = 1 and the identity function fi (x) = x. 
If there exists a Bernstein operator Bn fixing /o and fi, then the following equalities hold 
for the coefficients (3k and the nodes t^^k, whenever k = 0, n: 

(31) Pk = Pn-k and t^^k + tn,n-k = a + b. 

Proof. Let [/„ := {/(a + 6 — x) : / G Un}- Given g G Un, define Bng{x) := Bng{a + b — x). 
Since by symmetry Un = Un, and both operators i?„ and Bn fix the affine functions 
on [a,b], by uniqueness of the Bernstein operator we have Bn = Bn- Thus, Bnf{x) = 
Bnf{a+b—x) for every f E Un- Likewise, let pn,k, k = 0, ....,nhe a. non- negative Bernstein 
basis of Un- By suitably rescaling we may impose the following normalization: 



for all k = 0, n. Counting zeros at the endpoints, and using uniqueness of the Bernstein 
basis (up to a normalizing constant), we conclude that 















to ib) < ti ib) < t2 (b) < tg (b) =bforbe (0, tt) 
to {b) < ti (6) = 7r/2 = ta [b) < tg (&) = b for b = n 
to (b) < t2 (b) < ti (6) < t3 (b) =bforbe (tt, po) . 



(32) 




(33) 



'Pn,k (a + b-x) = Pn,n~k 



And now (IHT]) follows from (155]) and the fact that Bnf{x) = Bnf{a + b — x) for every 
feUn- □ 
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Proposition 27. Let pn,k, k = 0, ...,n, be a Bernstein basis, and for f G Un let (3q, ...,(3n 
be the coefficients in the expression 

n 

(34) / = 5^/3fcPn,fc. 

fc=0 

Then pn,n (b) = f ip) and 

(35) /' (6) = (3n-ip'n,n-l iP) + PnP'n,n iP) • 

Proof. For the first statement insert a; = 6 in (!34l) . for tlie second take tlie derivative of / 
in fl34l) and then insert x = b. □ 



Now we turn to the proof of Theorem | 
Proof. We define a Bernstein basis for U3 over [0, b] by setting 
P3,3 (a;) = X — sinx, 

P3,2 (x) = {b — sin 6) (1 — cosx) — (1 — cos 6) (x — sinx) , 
P3,i (x) = p3,2 {b - x) and ^3,0 (x) = ^3,3 (6 - x) . 

Let us check that this is indeed a Bernstein basis. We claim that the four given functions 
are positive on (0, tt) and have the required number of zeros at the endpoints. This is 
clear for p3 3, and hence for p^ Q. Regarding ^3 2, note that it has a zero of order 2 at 
and another zero at b. Since (1, x, cosx, sinx) is an extended Chebyshev space whenever 
b G (0, 2tt), P3 2 has at most 3 zeros in [0, b]. So it has no zeros in (0, b) and the zero at b 
is of order exactly one. Since P3 2(0) > 0, ^3^2 > on [0, b]. Thus, p^^i > on [0, b]. 

Recall that Bernstein bases are unique up to multiplicative constants. So to prove that 
a Bernstein operator does not exist, it is sufficient to consider the preceding basis. On 
the other hand, to prove that a Bernstein operator does exist, we need to exhibit nodes 
tk in [0,6] and positive coefficients a^, for k = 0,1,2,3. Now let 1 = ^^=0 '^'=^'3.^ ^^"^ 
X = J2k=o'ykP3,k- Lemma [5] tells us what the nodes and coefficients must be if Bs exists. 
By Proposition [27] we have 

Ps = 7 — and 73 ^ 



b — sin b b — sin b 

so ^3(6) := 73//33 = b. It follows from (15^ that 

(36) (32p's,2 (b) = -/53P3,3 (b) and 72^3^2 (^) = 1 - 73^3,3 (b) ■ 

Thus 

(b) = ^ = , 1 _ , _ b-sinb 

/32 -/33P's,Ab) -/33P'3,-Ab) l-cosfe- 

Now ( l3Til implies that 

b — sin b 
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We see that t2 (b) — ti (6) > if and only if g{b) := 2 sin 6 — 6 cos 6 — 6 > 0. Elementary 
calculus shows that g > on (0, tt), qItt) = 0, and g < {at least) on (tt, 3tt/2). If 6 = tt, 
then ti (tt) = t2 (tt) = tt/2. Furthermore, ^2 (b) < whenever sin 6 — 6 cos 6 < 0, so by 
Lemma [5l for 6 G {po, 27i) there does not exists a Bernstein operator. To see that such 
operator exists when b G (0,po); note that since /o = 1, by (fTT!) we have ak = Pk, so it is 
enough to show that Pk > for k = 0, 1, 2, 3. Since Po = /?3 and Pi = /?2 by Proposition 
l26l and /^s > 0, it suffices to prove that P2 > 0. Now from equation (136|) we get 



so > if and only if 6 sin 6 — 2 + 2 cos 6 < 0. Elementary calculus shows that this is 
the case for every b G (0, 27r), and in particular, for every b G (0, po]- 

Regarding the convexity assertions, if 6 G (0,7?] then the Bernstein operator preserves 
convexity by Theorem [22l Next, fix 6 G (7r,po]; write ti = ti (6), ^2 = ^2 {b)-, and consider 
the convex function / (x) = {x — ti) {x — ^2)- Since ti + ^2 = b, we have f {0) = f (6) = 
ti {b — ti) . By Proposition [26| /5o = /^s, so 



Using /3o/ (0) > we see that B^f is convex if and only if F := p^Q + 3 is convex. A 
direct computation shows that F (x) = b — sin {b — x) — sin x, so 

-F" (x) = sin (6 — x) + sin x. 

Thus -F" (0) = sin 6 < 0, since b G (vr, 27r) . By continuity, F" (x) < for all x in a small 



We see that non- increasing nodes can lead to pathological behavior on the part of the 
generalized Bernstein operator defined by them. Thus, either additional conditions are 
imposed to avoid this situation (in this case, a smaller value of b) or else the requirement 
that the operator fix 1 and x must be given up. On the other hand, there are natural 
examples where nodes are increasing, but not strictly. In addition to PropositionfTTj where 
the node at zero is repeated j times, we mention the case 6 = vr in Theorem [251 There, 
the second and third nodes are equal (to 7r/2). It is shown in [3] that the generalized 
Bernstein operator on U3 approximates some convex functions (such as |x — 7r/2|) on [0, vr] 
better than the standard polynomial Bernstein operator on (l,x, x^,x^) (while the latter 
operator approximates some other functions better). 
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neighborhood of 0, so F is not convex. 
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